Evidence for Thermally Activated Spontaneous Fluxoid Formation in 
Superconducting Thin-Film Rings 
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We have observed spontaneous fluxoid generation in thin-film rings of the amorphous supercon- 
ductor MoaSi, cooled through the normal-superconducting transition, as a function of quench rate 
and externally applied magnetic field, using a variable sample temperature scanning SQUID mi- 
croscope. Our results can be explained using a model of freezout of thermally activated fluxoids, 
mediated by the transport of bulk vortices across the ring walls. This mechanism is complementary 
to a mechanism proposed by Kibble and Zurek, which only relies on causality to produce a freezout 
of order parameter fluctuations. 



A number of years ago Zurek Q] proposed that con- 
densed matter systems can be used as models for test- 
ing ideas, first introduced by Kibble0,0|, concerning the 
quenching of topological defects in the early evolution of 
the universe. The principal picture developed by Kibble 
and Zurek was that, as a system approaches the critical 
temperature of a phase transition from above, the system 
slows down until at some temperature different regions 
cannot communicate with each other quickly enough to 
have correlated order parameters. As the transition pro- 
ceeds further, fluctuations between different regions are 
"frozen in", leading to relatively long-lived topological 
defects. Although this was only one of several mecha- 
nisms for freezout of topological defects considered by 
Kibble and Zurek, for brevity we will refer to this as 
the "Kibble- Zurek" mechanism in this paper. The ad- 
vantage of using condensed matter systems to test these 
ideas is clear: while it is impossible to do experiments on 
a cosmological scale, there are a number of well under- 
stood phase transitions into macroscopic quantum states 
in condensed matter that are experimentally accessible. 

The results of a number of experimental tests of these 
ideas have been mixed: Although vortex generation was 
at first observed in 4 He upon rapid quenching into the 
superfiuid state 0], a refinement of the original experi- 
ment by the same group found a residual vortex den- 
sity two orders of magnitude smaller than the Kibble- 
Zurek prediction. Experiments on vortex formation upon 
rapid quenching in 3 He 0, Q found good agreement 
with the Kibble-Zurek predictions. A search for sponta- 
neous vortex generation in a bulk superconductor upon 
rapid quenching @] was unsuccessful, although sensitiv- 
ity of 1CP 3 of the Kibble-Zurek prediction was reported. 
However, tests on superconducting loops interrupted by 
Josephson weak links have found good agreement with 
the Kibble-Zurek predictions, both for high-T c [9| and 
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superconductors. 
In addition to the standard Kibble-Zurek mechanism, 



which only relies on the satisfaction of the principle of 
causality to "freeze in" topological defects, it is also pos- 
sible for thermally activated defects to appear and be 
trapped in the system during cooldown pj . Detailed cal- 
culations of trapping of thermally activated fluctuations 
have been made in the hot Abelian Higgs model0|, and 
for the case of superconducting rings [lj]. In this paper 
we report on measurements of spontaneous fluxoid for- 
mation in arrays of thin film superconducting rings with 
various ring dimensions and compositions. SQUID imag- 
ing of such arrays allowed the performance of a number of 
cooling experiments simultaneously. The film thicknesses 
and compositions were chosen to obtain large magnetic 
penetration depths and inductances, which produce fa- 
vorable conditions for thermally activated behavior. The 
experiments provide evidence for the thermally activated 
appearance of topological defects (in this case vortices 
and antivortices) in systems undergoing a sufficiently 
rapid phase transition into an ordered state. We argue 
that such processes should be considered as a complemen- 
tary mechanism in tests of the Kibble-Zurek predictions. 

Although we obtained qualitatively similar results on 
rings fabricated from the low-T c superconductor Nb, and 
the cuprate superconductor YBa2Cu307_5, we will re- 
port in this paper on thin-film (G?=50nm) rings of the 
amorphous superconductor Mo^Si[l5T|. which showed the 
largest spontaneous fluxoid generation, and on which 
the most detailed measurements were taken. The films 
were fabricated by rf-sputtering on sapphire substrates 
and photolithographically defined using ion-beam etch- 
ing. Results will be presented here on two sets of rings. 
One set consisted of 25 circular rings, with inside diam- 
eter 2a=40 /Ltm and outside diameter 26=80 /im, spaced 
by 160 /im in a 5x5 square array. The second set had 
144 rings with 2a=20 /im and 26=30 /im, spaced by 60 
/im in a 12x12 square array. The films had a resistive 
superconducting critical temperature T c =7.81 K, with a 
transition width (10%-90%) AT C =0.03 K wide. Magnetic 
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FIG. 1: Scanning SQUID microscope images of a 12x12 ar- 
ray of 20 fim inside diameter, 30 pm outside diameter thin 
film rings of MoaSi, cooled in zero field through the super- 
conducting T c at various rates, and imaged with an 8 /im 
square pickup loop. The full scale variation of the grey scale 
images corresponded to 0.25"I>o (a), 0.3$o (b), 0.25 <3>o (c), 
and 0.37$o (d) change in flux through the pickup loop. The 
principal axes of the square array are rotated by 12° clockwise 
with respect to the image frame. 

imaging measurements were made with a variable sam- 
ple temperature, scanning SQUID microscope^fj, with 
a square pickup loop 8 /xm on a side. The sample was 
mounted with a thin layer of GE varnish on a saphire 
substrate to which a silicon diode thermometer and non- 
inductive resistive heater were attached. Both sample 
and SQUID were inside an evacuated can surrounded by 
liquid 4 He. The thermal time constant of the sample 
mount was varied by changing the amount of 3 He ex- 
change gas between it and the 4 He bath. Fits to SQUID 
images of the rings with fluxoid number N=l at zero 
field gave a Pearl length [l7| A=2A|/d (\ L the London 
penetration depth) =7 /im at 4.4 K. 

Our experiments consisted of cooling the sample in a 
given field at a controlled rate, and then imaging the 
rings with the SQUID microscope to determine the fi- 
nal fluxoid state. The spacing between rings was large 
enough that ring-ring interactions were not important in 
the cooldown dynamics. 

Figure 1 shows SQUID microscope images of the 20 
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FIG. 2: Probability of 20 ^m inside diameter, 30 /im out- 
side diameter MoaSi rings (closed symbols) and 40 /im/80 
fim rings (open symbols) having a final fluxoid number of +1, 
as a function of the cooling rate, after cooling in zero exter- 
nally applied magnetic flux. The dashed line is the prediction 
of the Kibble-Zurek quenched disorder model, and the solid 
lines are fits to the frozen thermally activated bulk vortex 
model discussed in the text. 

/im/30 /im ring array, cooled through T c at zero field with 
varying cooling rates. In this image the rings with final 
fluxoid number iV=0 are not visible; rings with positive 
fluxoid number (counter-clockwise circulating supercur- 
rents) are white; and those with negative fluxoid number 
(clockwise supercurrents) are black. Most of the rings 
with non-zero fluxoid number have N=+l or -1, but a 
few rings also have N=+2 or -2. For example, the upper 
left most ring in Fig. (d) has N=2. 

Figure 2 displays the results, plotted as the probability 
for a ring to be in the N=+l fluxoid state, for a number 
of cooling rates, for both ring sizes. The error bars are as- 
signed by setting the standard deviation of the number of 
N=+l rings to be the square root of the number of N=+l 
rings observed in each experiment. This plot shows that 
P(N=+1) varies roughly linearly with the logarithm of 
the cooling rate. 

The Kibble-Zurek prediction for the expected fluxoid 
number for a quenched superconducting ring due to free- 
zout of order parameter fluctuations is |l| 

N = ^VC/To(r/r Q )K (1) 

where C is the circumference of the ring, £o is the zero 
temperature coherence length, t = irh/lGksTc is the 
Ginzburg-Landau relaxation time, and Tq 1 =| dT/dt \ 
/T c is the normalized cooling rate. The dashed curve 
in Fig. is this prediction for the 20 /xm/30 fim rings, 
using £o=4-9nm 18], and taking C to be the inside cir- 
cumference of the rings. The corresponding prediction 
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for the 40 /mi/80 /im rings would be scaled vertically by 
a factor of y/2. The Kibble- Zurek mechanism predicts 
more spontaneous fluxoid generation than is observed, as 
well as a strong increase with ring size which is not ob- 
served. Moreover, the Kibble-Zurek prediction, although 
it has a weak (Tj 1 IdT/dtl) 1 / 8 dependence, nevertheless 
varies more rapidly with cooling rate than experiment. 
We should note that the Kibble-Zurek predictions for a 
superconducting ring were for ring walls much thinner 
than the ring diameter 0. However, we do not believe 
that our finite ring wall widths would effect the Kibble- 
Zurek predictions qualitatively. 

We suggest that our observed logarithmic dependence 
of non-zero fluxoid probability on cooling rate is due to 
spontaneous changes in the ring fluxoid number by ther- 
mal activation in, and transport across, the ring wall 
of bulk vortices. Such a model has been used to qual- 
itatively explain measurements on telegraph noise and 
fluxoid transition rates in superconducting thin-film rings 
of underdoped cuprate BioSroCaCuoOe_s[lfll.l2(Tj. which 
have penetration depths and coherence lengths compara- 
ble to the present samples. Physically our picture is that, 
as the rings cool through T c thermally activated fluctu- 
ations rapidly change the fluxoid number. The thermal 
activation barrier to this process becomes large as the 
temperature is reduced further, eventually turning this 
process off, leaving the rings with some probability of 
being in energetically excited states. 

We calculate the probability Pn of a ring being in the 
Nth fluxoid number state as follows: Consider a set of 
ring fluxoid states N, with ring energies En, and with 
barrier energies Eb (in our model the energy required to 
nucleate a vortex in the ring wall) that must be overcome 
to make the transition between them. The transition rate 
v n ^n> is given by 



P N ^e- {EB - EN)/ki 



(2) 



where v$ is an attempt frequency, and we neglect the de- 
pendence of the transition rate on the final state energy. 
At high temperatures, in equilibrium, vn^n' = vn'^n 
for all N, N'. This implies that 



Pn =e 



E N /k B T P -E N , /k B T 



N' 



(3) 



The ring energies E N = E r (N - <f> a ) 2 , with <f> a = $ a /$ 
the normalized externally applied flux, and $o = hc/2e 
the superconducting flux quantum[2]| . Eq. 03 implies 
that the probability of a ring having a non-zero flux- 
oid number depends on the ratio E r /k B T. As the 
ring cools down, the transition rates slow down when 
Eb(T) ~ fcflT, fixing the ring fluxoid number. We can 
assign a "freezing" temperature T> by insisting that the 
probability of there being a further fluxoid jump, ob- 
tained by integrating the transition rates over tempera- 
ture from T=Tp to T=0, be less than 1. It can then be 
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FIG. 3: Probability of 20 /im inside diameter, 30 /im outside 
diameter M02S1 rings having a final fluxoid number of +1, as 
a function of the externally applied magnetic flux. The open 
symbols are the data, the solid line is a fit to the thermally 
activated vortex model discussed in the text. 



shown that 

k B T F 



-E B {T F )/hx 



7 I dT/dt 
k B v§Tp 



(4) 



where 7 = dE B (T)/dT \ t =t f - We then rewrite Eq. El to 
obtain the probability Pn,f that a ring is frozen into a 
state with fluxoid number N: 



P. 



N,F 



7£' 



-Xr(N'-<P a ) 2 



(5) 
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with Xr = E r (Tp)/k B Tp. For the parameters of our 
system, Pn,f varies nearly linearly with \r, an d therefore 
depends roughly logarithmically on cooling rate through 
Eq. H 

The upper solid line in Fig. [S] is a fit of Eq. [5] to 

the 20 /im/30 /im ring data, resulting in best fit values of 
£ r (T F )/.E B (T F )=0.097±0.1and 7//c B ^ =2.4x 10~ 7 sec, 
with upper and lower limits of 17.5xl0~ 7 sec. and 
0.4xl0 -7 sec. respectively on the latter value, using 
a doubling of \ 2 as a criterion. In the temperature 



range of interest 
A(T) ~ A(0)/(1 ^ 
sions. In this limit 



1-t: 

4' 



imit|2q| 



=(T C -T)/T C < 1, the Pearl length 
is larger than typical ring dimen- 



E r {A > a,b) ~ $g/(87r 2 A). 



(6) 



The energy E B required to nucleate a vortex in a thin 
film strip of width W is given by (2^ 



E B = 



ln(— 

8tt 2 A n V tt£ 



(7) 
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Combining Eq.'s HO and we estimate E r /Es ~ 
l/ln(2H70. Usin g Eq. 0, Eq. Q, A(0)=7 /jm, 
£(T) = ^(O)/^! - 1, and TV = 6 - a = 5 /mi, leads 
to (T c - T F )/T C = 2xl0 -3 , 7 =2.7xl0 -14 erg/K, lead- 
ing to EV/-E.B =0.30, a factor of three larger than our 
fit value. Our fit value for implies t^o=8.3xl0 8 

sec." 1 , with upper and low limits of 6.1xl0 9 sec. -1 
and 1.1 x 10 s sec. -1 respectively, in reasonable agreement 
with attempt frequencies obtained from fluxoid dynam- 
ics of Bi^SroCaCu^Os-A rings with similar penetration 
depths Qjltij. The lower solid line in Fig. H is a 
fit to the 40 /zm/80 /tm ring data, using the value for 
j/kgVQ obtained from the 20 /tm/30 /tm rings, and al- 
lowing E r (Tp) I Eg(Tp) to vary. The best fit value in this 
case is E r (Tp) / Eb(Tp)=0.12, with upper and lower lim- 
its of 0.15 and 0.105, to be compared with a calculated 
estimate for these rings of E r /EB=0.2. Note, however, 
that the predicted dependence on wall width does not 
match the experimental results: doubling the size of the 
ring wall should decrease the value for E r /Es by 30%, 
but our fit values increase by about that amount. 

A further test of our model is provided by measure- 
ments of the final state fluxoid number probability as 
a function of applied field. Figure 01 shows the results 
of such experiments for the 20 /tm/30 /tm rings, cooled 
at a cooling rate of 15 K/sec. At first these results are 
startling: the fluxoid number probabilities depend only 
weakly on applied field. This is in contrast to pioneering 
work by Davidovic et al. |23j,|24J on 1.6 /tm outside diam- 
eter thin film aluminum rings, which showed, e.g. Pi.f 
stepping up sharply from to 1 at </> a =0.5. This differ- 
ence can be explained qualitatively as follows: our rings, 
which have large values for a, b, and A, have small ring 
energy spacings (scaled by E r ). Therefore ksTp is com- 
parable to E r (Tp) and it is probable that energetically 
unfavorable final fluxoid numbers will be populated. On 
the other hand, the rings of Davidovic have much smaller 
values for a, b, and A, and the ring energy spacings will 
be correspondingly larger. The solid line in Fig. [21 is a fit 
to the data using Eq. |SJ with \r as a fitting parameter. 
The best fit is obtained for \ r = 1.75, compared to a 
value of 4.3 calculated from Eq. 0] 

In conclusion, we have found that our results for the 
spontaneous generation of energetically unfavorable final 
fluxoid numbers in cooled thin-film rings of the amor- 
phous superconductor MoaSi can be explained using a 
model of thermally activated bulk vortices. This mecha- 
nism occurs over a much different time and temperature 
scale than the conventional Kibble-Zurek mechanism, 
which predicts that freezout of order parameter fluctu- 
ations occurs at e =| (T-T c )/T c \= ^t /t q wlx 10 -7 . 
For our rings the thermally activated processes persist to 
e w 1 x 10 -3 . In this case any fluxoids generated by the 
conventional Kibble-Zurek mechanism would be "washed 
out" by the thermally activated vortex mechanism. Our 



rings are exceptionally susceptible to thermally activated 
vortex processes: they have very large magnetic penetra- 
tion depths and inductances, which make nucleation of 
vortices not too energetically costly at temperatures rel- 
atively far from T c . Nevertheless, the possibility of ther- 
mally activated processes should be considered, keeping 
in mind that not only short coherence lengths, but also 
short penetration depths, are desirable when using su- 
perconductors for tests of the Kibble-Zurek mechanism 
for topological defect formation in quenched macroscopic 
quantum systems. 
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